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We investigate the role of noise in the phenomenon of stochastic synchronization of switching 
events in a rocked, overdamped bistable potential driven by white Gaussian noise, the archetype 
description of Stochastic Resonance. We present a new approach to the stochastic counting process 
of noise-induced switching events: starting from the Markovian dynamics of the nonstationary, 
continuous particle dynamics one finds upon contraction onto two states a non-Markovian renewal 
dynamics. The output frequency is determined as the velocity of the underlying discrete phase 
dynamics. The phenomenon of noise-assisted phase synchronization is investigated in terms of an 
effective, instantaneous phase diffusion. The theory is applied to rectangular-shaped rocking signals 
O ' versus increasing input-noise strengths. Precise numerical simulations corroborate very favorably 

our analytical results. The novel theoretical findings are also compared with prior findings. 
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^ . I. INTRODUCTION 

ts. 

The theme of synchronization has wide-spread applications, covering a plenitude of phenomena 0, 0, 0, 0] • Some 
characteristic examples are the entrainment of a system by external, time-dependent forcing, or thegeneralization of 
the synchronization concept for systems that exhibit a chaotic dynamics [j| , lag synchronization U and also phase 
synchronization 0. Synchronization phenomena play not only a key role for diverse technological applications but 
increasingly as well for the description, the control, and even for the therapy of selected medical disorders @ . 

Due to the interaction with a surrounding environment or with internal degrees of freedom, noise is present in many 
physical systems. This being so, its role cannot be ignored when investigating synchronization phenomena. In recent 
years, it has turned out that noise can actually play a constructive role in many physical situations. In particular, 
noise can boost the transduction of information by means of the phenomenon of Stochastic Resonance [9j in an ample 
number of metastable physical and biological systems [To[ . Furthermore noise enables Brownian motors to do work 
against external load forces , or to induce phase transitions far away from thermal equilibrium [l2L Il3| . 

Our focus here is on the role of phase synchronization in stochastic overdamped systems driven by white Gaussian 
t— I ' noise. In these cases, the velocity of the dynamics is not a measurable quantity because the stochastic trajectories 
are neither differentiable nor of finite variation, see e.g. Ref. 14]. A recently proposed method for measuring the 
average phase velocity or frequency which is based on the generalization of a Rice rate formula for threshold crossings 
is consequently not a suitable method [l5|. An alternative approach is based on the so-called " Hilbert phase" 
dynamics, as pioneered by Gabor [Tfjl ] for deterministic systems. In the present work, we shall take a closer look at the 
synchronization phenomenon in a periodically driven bistable system. Then it is advanta geo us to introduce a discrete 
phase dynamics, as recently proposed by Schimansky-Geier and collaborators 0,0,0,13- In order to extract this 
discrete phase dynamics from the underlying continuous process, we shall consider the stochastic counting process 
of the noise-induced switches between the two potential minima. It turns out that this counting process is in fact a 
O ■ non-stationary renewal process 0] . 

The outline of the present work is as follows: First, we introduce a dichotomic process by filtering out the fluctuations 
around the potential minima of the original stochastic process, which possesses a clear interpretation in terms of 
a discrete phase. By contrast to the underlying stochastic process, this two-state process, however, is no longer 
Markovian. Subsequently, in Sect. Ill, we analyze in detail the statistical properties of the random switching-times 
associated to the dichotomic process. The one-time statistical properties of the discrete phase are then studied in 
Sect. IV. Based on these results, exact analytical expression for the instantaneous output frequency and the phase 
diffusion are derived. Approximate expressions are then derived for the weak-noise limit and for a slow external 
driving. Finally, our analytical findings will be applied to the case of a symmetric bistable potential driven by a 
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periodic rectangular input signal. To corroborate our analytical results, we compare them with those obtained from 
a numerical simulation of the original stochastic process. 



II. DESCRIPTION OF THE MODEL AND DEFINITION OF THE DISCRETE PHASE 

To start, we consider a stochastic dynamics characterized by a single degree of freedom x(t), whose dynamics (in 
dimensionless units) is described by the stochastic differential equation 

x(t) = -U'[x(t),t]+£(t), (1) 

where is a Gaussian white noise of zero mean with autocorrelation (£(i)£(s)) = 2D8(t — s), and U'(x,t) is the 
derivative with respect to x of the bistable quartic potential 

U(x,t) = ^-^-F(t)x, (2) 

F(t) representing a periodic forcing with period T. Our focus is on subthreshold signals; more precisely, we will 
assume that, for any instant of time, \F(t)\ < A^ — 2/V27 where A t h is the static threshold value (the dynamical 
threshold value always exceeds this adiabatic threshold Ath)- In this case, the potential possesses two minima at 
q~i(t) < and q+i(t) > 0, and a maximum at qM (t)- Introducing the function v(t) = arccos [F(t)/ A t h], with arccosy 
being the principal value of the arc cosine of y (i.e., the value in the interval [0, 7r]), and 
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= 71 cos 



v(t) + 2wn 



(3) 



then rjo(t) yields the location of the minimum to the right of the barrier [i.e., q+i(t) — T)o(t)], rji(t) yields the location 
of the minimum to the left of the barrier [i.e., q~i{t) — r)i(t)], and r]2(t) yields the location of the maximum [i.e., 
qAi(t) = %(£)]■ From now on, we will assume that at an initial instant of time to the system is placed at one of 
the minima of the potential q ao (to), with cto = +1 or — 1. The long-time behavior of the quantities of interest can 
be obtained by taking the limit as t — > — oo at the end of the calculations. Henceforth, we will make explicit the 
dependence of all the quantities on the initial preparation by the superscript ao>*0) see also discussion after Eq. il4*0")l . 
Thus, for instance, we will write x a °' ta (t) instead of x(t), meaning that x a °' ta (to) = q ao (to). 

To analyze the synchronization phenomenon in this stochastic bistable system, it is convenient to introduce a 
discrete phase associated to the continuous stochastic process x ao,to (t). In order to do so, firstly we will proceed 
to filter out the fluctuations around the minima of the stochastic process x ao,t °(t) to obtain a two-state stochastic 
process x a °' to (t) which only takes the values +1 or — 1. The procedure used is as follows: At the initial instant of 
time to we set x a °'*°(*o) = c*o- A switch of state from ±a to =F«o occurs whenever the system, having started in one 
of the minima, reaches the other minimum for the first time. The instant of time at which the n-th switch of state 
takes place is a random variable which will be denoted by T" 0,t °, with n= 1,2.... Formally, these random variables 
can be defined recursively as 

T„ Qo < to = min [t : t > T*"? and x a °> to (t) = q an (i)] , (4) 

where T Q a °' to = t and a n = (-l) n a . Thus, if we introduce the stochastic process 

N ao - to (t) = max [n : T n Q °'*° < t] , (5) 

which counts the number of switches of state in the interval (to,t], then the two-state stochastic process x a °''°(i) can 
be expressed as 

x a °' to (t) = a cos [irN a °' to (t)] . (6) 

By analogy with the case of a sinusoidal signal, we will define the discrete phase ip a °' to (t) associated with x a °' to (t) as 
the stochastic process 

ip ao > to {t) = TrN ao - to (t). (7) 

In Fig. ^ we illustrate the procedure just described in the particular case of a rectangular signal [see Eq. I|tj5fl and 
text below] with amplitude A = 0.25 and frequency u> — 2tt/T = 0.01, and a noise strength D = 0.02. The initial 
instant of time has been chosen to be to = 0, and the system has been initially placed at q+i(0), so that ao = +1. In 
the upper panel, we have sketched the rectangular periodic signal F(t), whereas in the middle and lower panels, we 
have depicted a random trajectory of the stochastic process x +1 '°(t) and the corresponding realization of the process 
respectively. 
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FIG. 1: Illustration of the procedure used to define the stochastic process x" '' W i n the particular case of a rectangular signal 
[see Eq. 1651 and text below] with amplitude A — 0.25 and frequency SI = 2n/T — 0.01, and a noise strength D — 0.02. The 
initial instant of time has been chosen to be to = 0, and the system has been initially placed at 5+1(0), so that ao = +1. In the 
upper panel, we have sketched the rectangular periodic signal F(t), whereas in the middle and lower panels, we have depicted 
a random trajectory of the stochastic process x +1 '°(t) and the corresponding realization of the process X +1 '°(t), respectively. 



III. STATISTICAL CHARACTERIZATION OF THE SWITCHING TIMES T"°' t0 

According to Eqs. JSJ and J7J), the statistical properties of the discrete phase ip a °' to (t) are closely related to those 
of the random switching times T"°' to . The aim of this section is to provide a detailed description of the statistical 
characterization of these random variables previous to the analysis of ip a °' to (t), which will be postponed to the next 
section. The connection between the statistical properties of the switching times and the original stochastic process 
x ao,t °(t) will be also analyzed in this section. 

Following the approach presented in Ref. [l^, the random variable T^ 0it ° can be characterized statistically by its 
probability density function 

tf^W- A M " < t^" 1 . (8) 

At^0+ At 

Besides this probability distribution function, it is also convenient to introduce the (cumulative) distribution function 

G%"*°(t) = Prob [T n a °' io <t] = f dt' g%°' to (t'), (9) 

as well as its complementary, 

g%°' to (t) = Prob [T„ Qo < to > f] = 1 - G%°' to (t). (10) 

In the particular case n = 1, these functions can be directly determined from the solution of the Fokkcr-Planck 
equation (FPE) 



-P(x t) = — 

dt [ ' ' dx 



D^- + U'(x,t) 



P(x,t), (11) 



with initial condition 

P(x ) t ) = 6[x-q ao (t )}, (12) 
and absorbing boundary condition at q^ aa (t), i.e., 

P[q^ ao (t),t] = for alU > t . (13) 
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Denoting by P Q °>*° [x, t) the solution of the above problem, it follows from the definition of Gi°' to {t) in Eq. (|10|) . with 
n=l, that 

g?°> t0 (t)=6 ol0t+1 dxP+^ix^+S^x dxP- l ' t0 {x,t), (14) 

J q—l (t) J — oo 

for t >to, and = 1 for t < to. The function ' (t) is the conditional survival probability of the discrete state 

«o- The function G"°' to (t) is then given by Eq. (|10fl with n = 1, whereas gi° ,to {t) = — Gi°' °(t) is the corresponding 
conditional residence time distribution (RTD). The dot indicates the derivative with respect to time t. The knowledge 
of either the conditional survival probabilities G""' °(t), or (equivalently) the conditional RTDs g" ' °(t) is sufficient 
to specify a driven two state non-Markovian renewal process |20| . These functions can be found from the underlying 
continuous state Markovian dynamics by solving Eqs. Ijll|l - (|14|l . 

For n > 1, the functions &" '*°(i), G"°' to (t), and can be obtained iteratively from the ones corresponding 

to the case n = 1 by making use of three integral equations. To obtain the first integral equation, let us consider the 
consistency condition 

/+oo 
dt' Prob [T«°f° >t\%?<" to = t'] gT M (t'). (15) 
-oo 

By use of the definition of the switching times in Eq. (0J , as well as the Markovian character of the original stochastic 
process x a "' to (t), it is straightforward to verify that 

Prob [T£°'*° > t 1 T,^ t0 = t'] = Gi^' (t). (16) 

Inserting the above expression into Eq. I|15|l and taking into account that (t) = 1 for t < t' and S , ™ 0,t °(i) = for 

t < tg, we obtain 

Ci°(t) = G?'**® + f dt' G^'^g^it'), (17) 

Jto 

for n > 1. The interpretation of this result is straightforward: The probability that the (n + l)-th switch of state 
occurs after the time t is equal to the probability that the n-th switch occurs after that instant of time plus the 
probability that the n-th switch has happened at any time t' before t with the next switch taking place after t. 
Similar interpretations hold for the integral equations; i.e., 

GTA° (*)= f dt' G^'WC^i'), (18) 



Cx°(*)= / dt' g^'it^it'), (19) 



which are obtained from Eq. I|17|) by using Eq. (|10fl and g"°' to (£) = — Gn 0:t °(t), respectively. 

A formal solution for g" 0, *°(i n ) is obtained by solving iteratively the integral equation (|TT)|) . The result is 

9n°' ta {tn)= dt n . x ... dh Hg^(t j+1 ), (20) 

J t Jt j = Q 

for n > 2. Thus, the probability distribution corresponding to the first switch of state, <?" '* (i), which can be obtained 
from the solution of the FPE (|llfl with i|12fl and i|13fl , determines completely the statistical properties of the rest of 
the switching times. 

IV. ONE-TIME STATISTICAL PROPERTIES OF THE DISCRETE PHASE: THE OUTPUT 
FREQUENCY AND THE PHASE DISPERSION 

The one-time statistical properties of the discrete phase ip ao,t °(t) can be evaluated by making use of the probability 
distribution of the number of switches of state 

p%°' to (t) = Prob [N ao < to (t) = n] , (21) 
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with n = 0, 1,2, .. .. From the definition of N a ^ to (t) in Eq. 10 it follows that 

Prob [N ao ' to (t) >n]= G%°> to (t), (22) 

with Go '*°(£) = 1. Consequently, the probability distribution of the number of switches of state and its derivative 
with respect to t can be expressed, respectively, as 

PT tQ (t) = G?*(t)-G°$>(t) 

= G: o +i o (t)-GT to (t), (23) 

and 

Pn°' t0 (t)=g^ t0 (t)-g:° + t 1 °(t), (24) 

with £o 0,to (i) = 9o ' to (t) = 0. The average of an arbitrary one-time function of N a °' to (t), K [N a °' to (t)], is obviously 
given by 

oo 

(K [N<*»*°(t)]) =Y,K{n)p^°{t). (25) 

Equations (|23|l and (|24|l can be written in a more transparent form by introducing the probability of an almost 
immediate switch of state after n switches 

Prob It <T"°i*° < t + At I N ao ' to (t) =n] , x 

r"°< to t = lim L - n+1 ~ 1 ^ 1. 26 

n K ' At^o+ At v ; 

Note that if the process x a °' to (t) were Markovian, these probabilities could only depend on the state a n . The explicit 
dependence on the number of jumps n and on the initial preparation ao at time to is a consequence of the non- 
Markovian character of the process. Another fingerprint of the non-Markovian nature of the dichotomic process is the 
fact that these probabilities depend on the time t even in the absence of the external driving. In order to clarify this 
point, let us consider, e.g., the particular case n = 0. Then, while initially, right after the particle has been prepared 
at one of the minima, the distribution function P a °'*°(a;,t) is still very sharply peaked around q ao (to), it becomes 
smeared out around the minimum after the intra-well relaxation time. Consequently, the probability of an immediate 
switch is obviously will be different in both situations, even without an external driving. 

Multiplying and dividing the right-hand side of the above expression by p" 0: *°(t) and taking into account that 



Prob [t < T"^° < t + At and N a °- to (t) = n] 
At^o+ At 



linr 1 n + i - _ w L = c m {t) (27) 



it is readily seen that 

Then, it follows from Eqs. (JTJl, (EH, and (gBJ that 

p^°(t) = fdt 1 sr"<*'(t)r^f (t>^f (t'), (29) 

Jt 

for n > 1, with Po°' to (t) = Gi°' to (t)- Analogously, Eq. (|24|l leads to the following hierarchy of differential equations 

Pn°' t0 (t) = K-l(t)pT-l (t) ~ K > t0 (t)pT t0 (t), (30) 

for n > 1, and 

Ar io (*) = -C' to (%o°' to (*), (3i) 

which must be solved with the initial condition /?" '*°(to) = 5 n ,o- 
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We will also introduce the conditional probability for x Q °'*°(0 to take the value j3 = ±1 at time t, provided that it 
took the value ag with probability one at the initial instant of time to, 

^• t0 (t) = Prob[x Q °' t0 (*)=/?]. (32) 

Noting that after an even number of switches of state the system ends up in the same state as it was initially, whereas 
for an odd number of switches the system ends up in the other state, it is clear that the events {x a °' to (t) = a o} 
and {x a °' to {t) = — «o} are respectively equivalent to the events {N ao,t °(t) is even} and {N a °' to (t) is odd}, and con- 
sequently 

oo oo 

P a p°' t0 it) = S ao , £ P 2 M '(t) + 5. ao ,p pZ+M (33) 

n=0 n=0 

Besides the probability distribution of the number of switches of state, later we will also use the probability 
distribution of the number of switches of state conditioned to the value of X a "' to (t)i 

p%°' to (t\f3) = Prob [N ao ' to (t) = n | X ao,t °(t) = 0\ . (34) 

Multiplying and dividing the right-hand side of the above expression by p'g ' °{t) and taking into account that 
Prob[A Q °<*°(i) = n and x"°' to (£) = P] = p" ' t °(0 <S a„, / 3, it results that 

^°' t °(^) = ^|i||^, /3 . (35) 
Pp \ l ) 

The average of an arbitrary one-time function of N a °' to (t), K [N a °' to (t)], conditioned to the event {x ao,t °(t) = (3} 
will be denoted by (K [N a °' to (t)]) /3 , with (3 = +1 or - 1, and it is given by 

OO 

n=0 

oo a a ,t /.n, oo a ,t (,\ 

= <W E K <M%*M + '-o* E + (36) 

n=0 P<*o W n=0 P-ao W 

^From the above expression and Eq. 1(25(1. it follows that 

(A" [N ao ' to (*)] ) = (K [N ao ' to (*)] > +1 plf° (t) + (K [N ao ' to (t)] )_ x p a _f° (t). (37) 

Another interesting quantity which will be useful later is the probability of an almost immediate switch from state 
(3, defined as 



Prob 



t<T%°< % {t)+1 <t + At\ X a °' t °(t)=(3 



i?«(t) = A hm + L • (38) 

Notice that in the Markovian limit these probability are independent of the initial preparation ao at time to and 
cannot be distinguished from the probabilities defined in Eq. 1)26(1 . Multiplying and dividing the right-hand side of 
the above expression by pZ ' °(t) and taking into account the equivalence of the event {x a °' to (t) — (3} with the event 
{N a °' to (t) is even} 5 Qo>/ 9 + {N a °> to (t) is odd} 8- a0tf3 , as well as Eqs. (EU and it is easy to see that 

7|°' t0 (*) = (r^4%,(t)(*))^ (39) 

Differentiating Eq. 1(33(1 with respect to t and taking into account Eqs. 1(30(1 . 1(36(1 . and 1(39(1 . it is straightforward to 
obtain that 

pT t0 {t) = -^°' to (t)p a e ' to (t) + ^f°(t) P a _f°(t), (40) 

for (3 = 1 and — 1. The Eq. 140(1 is a non-Markovian master equation for the conditional probabilities P^j°' to (t)- It 

is of the time-convolutionless form 0,H3. The rate parameters 7?°'*°(t) entering this equation are time-dependent 
quantities even in the absence of time-dependent driving. Such a time-dependence reflects primarily a non-exponential 
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distribution of the residence times of the renewal two state non-Markovian process |19l |23( which results from the 
projection of a continuous state Markovian stochastic dynamics onto the two discrete states j3 = ±1. A time-dependent 
driving introduces an additional time dependence into 7? ' °(t) which is present also in the driven Markovian case. In 

this latter case, 7 i g 0,t °(i) becomes a time-dependent rate and this rate depends neither on a®, nor on to (see below). 
Two other forms are possible to describe the evolution of conditional probabilities Pg°' to {t). One is given by the 
generalized master equations (GMEs) with the memory kernels expressed via the corresponding RTDs. In a driven 
case, the kernels of corresponding GMEs will become functionals of driving and will depend on both time arguments. 
Alternatively, the integral equations for the conditional probabilities p^j°' ta (t) can be derived for the driven two-state 

renewal process in terms of (conditional) RTDs <7i°' io (t) [2(J. Such integral equations present a generalization of the 
integral renewal equations of Ref. [191 onto the driven case. We apply in this work a time-convolutionless description 
of non-Markovian dynamics |2lL l22j| to the synchronization problem. 

After these rather formal considerations, we shall now apply these results to the evaluation of two important 
quantities in the study of the synchronization phenomenon: The instantaneous output frequency and phase diffusion. 



A. The instantaneous output frequency 

The instantaneous output frequency is defined as the derivative with respect to t of the averaged discrete phase 
[13, i.e., 

<Ct t0 (t) = I (^ t0 (t)) = tt| (N****®) . (41) 

Multiplying Eq. (EUJ) by 71, summing up the series X^^Li np" 0,t ° (t), and taking into account Eq. (|25|l . it is easy to 
obtain that 

§ i (N ao ' t °(t)) = (rX\ {t) (t))- (42) 

Consequently, from Eqs. Ip?7| . J2HJ), and lHT|) . it results 

fCt'° (t) = tt [73'*° (t) + 7 ^ t0 (t)plf (t)] . (43) 

This finding for the averaged frequency of the discrete phase dynamics constitutes a first main finding of this work. 



B. The instantaneous phase diffusion 

Let us now proceed to the evaluation of the instantaneous phase diffusion D"°l to (t), which can be defined as 0] 

d 



D 



= ^{([<° , *°( i )] 2 )-<¥' Q!0 ' t0 (*)> 2 } 



^|{(r°' fo w] 2 )-( ro,io w) 2 } 



(44) 



Multiplying Eq. by n 2 , summing up the series X^^i n2 Pn 0,t °(t)i anc ^ taking into account Eq. I|25|l . it is straight- 
forward to see that 

I ([N<*°>t°(t)] 2 ) = 2 (N^(t)T a X\ {t) (-t)) + (r#&o (t) (*)) • (45) 
Replacing the above expression into Eq. (|44|l and taking into account Eqs. I|41|l and (|42|l . it results that 



-2-k 2 



N a ^{t)T%°*\ {t) {t)) (N a ^{t)) (r^4° t0(t) (i) 



(46) 



Equation (|46|l can be expressed in a more convenient form by writing all the averages (...) in term of the conditional 
averages (. . . )p, according to Eq. (|37|l . Then, after some simplifications one obtains 



DZ?°{t) = TittZt a {t) + 27r 2 A7 Qo ' to (0* Qo ' to W 

+ 27T 2 Y, Cl°' t0 {t) P f^{t), 



(47) 



0=±1 
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where 



A7 Qo ' to (t)=7+r o (*)-7-r°(*), 



(48) 



* Q0 ' t0 (t)= {N a ^ fo (t)) +1 - (N ao ' to (t))_ 1 
and we have introduced the conditional covariance 



Cp oM, (t) = ( N a °- ta (t) - (N a °- ta (t)) 



N"o.to(t)\ L ) \ L N"o,to(t)\ L ) / 



(49) 



(50) 



Making use of Eqs. l|3Ufl and l|4(J|l and after some lengthy calculations, it is possible to prove that \fr a °' to (t) satisfies 
the differential equation 



ir ao ' to (t) = -j ao - to (t)y ao < to (t) 



0=±i 



where 



7 «»*(()= 7 ;f(i)i 7 r(i). 



(51) 



(52) 



Equation (|5f \ can be formally solved taking into account that, as it follows from the definition (|49|) . \I/ Qo ' to (to) = 0. 
The result is 



p=±i 



to 



x^°' to (t')e-/*' dt "^ 0,t °( t "). 



(53) 



Replacing the above expression into Eq. (|47|l . one obtains 
D«f°(t) = ^^ to (0 + Q Qo ' to W 



J t t / d*"7 Q 0' t0 (t") 



(54) 



where all the dependence on the conditional covariance C? ' °(t) has been included in the function 
Q ao > to (t) = 2tt 2 ^ C , ;°' t0 (^°' t0 (i) 



2 7r 2A 7 Q! °'*°(i) E ^ / ^'C; 0! * (i , )P^'*°(i')e" /tt " : ' t " 7Q0 ' t0(t " ) - 



(55) 



The expresion (|54|l for the instantaneous phase diffusion presents a second main result of this work. 



C. The weak- noise and low- frequency limit 

Throughout the following, we will assume that the noise strength D is sufficiently small so that the intrawell 
relaxation time scale is negligible compared with the time scale associated to the interwell transitions and, as well, 
the driving time scale T. In this case, for t — to much larger than the characteristic intrawell relaxation time, the 
probability of an almost immediate switch of state after switches, TQ' to (t), can be approximated by the Kramers 
rate of escape [24| from the state (3 at time t, i.e., 



ir°(i)«7*(*) = 



W/3 (t) UJ M (t) / U [q M (t),t]-U{qp(t), t] 



2tt 



■ exp 



D 



(56) 
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where oj p {t) = y/U"[q p (t),t] = ^/3[qp{t)} 2 - 1 and w M (t) = y/\U"[q M (t),t]\ = ^1 - 3[q M {t)} 2 . Furthermore, from 
Eqs. (|TT|) . (J2SJ, and (J2HJ) it follows that, within this approximation, we also have that r^'*°(t) s» 7g(i) for n > 1 

and, consequently, jg°' to {t) « 7^(i) and C^°' to (t) « 0. In this case, the conditional survival propabilities and the 
residence time distributions in state ao read 



(57) 



and 



<?r' t0 W=7^ Wexp - / 7 £(f)d* 



(58) 



respectively. This corresponds to a two-state Markovian process with rates 7±x(£). In this Markovian limit for the 
reduced, two-state dynamics, the instantaneous output frequency and phase diffusion become, 



*Ct*°(*) = * [7+1 Wf^W + l-i(t)P-°i t0 (t)} 



(59) 



and 



27r 2 A 7 K (0 ^ /} / dt'^{t')\p^{t') 



3 -/ t Ut" 7 K (t") 



(60) 



respectively. Here, 7 K (i ) = 7+1 (t) + 7^ (i) , A7 K (t) = 7^ (i) — 7?i (t) , and ,to (t) is obtained by solving the master 
equation 



pf<\t) = -i&{t)pf™{t) +7- /3 WP-°/°(*), 



„K/,\ o ,fo 



(61) 



with initial condition p^°'*°(to) = Sa ,f3- 

In order to obtain expressions independent of the initial preparation, it is necessary to take the limit to — ► —00 
of Eqs. I|59(l and (|60|l . In this limit, it can be shown that the functions f2 out (i) = lim^^-oo fl"°f°(t) and D out (t) = 
Imitg^-oo D °i °(t) are periodic functions of the time t. Then, one can perform a cycle average and define the averaged 
output frequency 



- I dtn out (t) = J / dt [7?i(t)p+i(t) +7-1WP-1W] 



r 



(62) 



wherep^(i) is the periodic long-time solution, pp (t) = lim to ^_ 00 p^ '*°(f), ofEq. (|B*T|l . After some lengthy calculations, 
it is also possible to show from Eq. (ffiTl)) that the averaged phase diffusion is given by 



T 



Amt = TP dtD out (t) = irn out - — csch 



T 



7 K T 



P [ T dt f ' dt'A^it^fWlppit'rfe^ 
0=±i Jo Ja 



(*-*') ^-/;*"7 K ('") 



(63) 



with 



7 K = ^/ *7 K W- 



(64) 



In the next section we will consider the case of a rectagular input signal. In this case, explicit analytical evaluations 
of the integrals in Eqs. I|62[l and (|63Jl can be carried out. 
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V. PERIODIC RECTANGULAR INPUT SIGNAL 



As an example of the use of Eqs. I|62[) and l|63() which is amenable to analytical treatment, we will consider the case 
of the periodic rectangular driving force 



F(t) = (-!)"« A, 



(65) 



where n{t) — [2t/T\, [z\ being the floor function of z, i.e., the greatest integer less than or equal to z. In other 
words, F(t) = A [F(t) = —A] if t £ [nT/2, (n+ l)T/2) with n even (odd). Because the potential fulfills the symmetry 
property U(x, t + T/2) = U(-x,t), we have q M (t) = (-l) n Wg M (0), and 

^ = (-iy^-(-ir^, (66) 

where Aq(Q) = q+i(0) — <7_i(0). Here, we have taken into account Vieta's formula q\{t) + q^t) + quit) = 0. 
According to the above mentioned symmetry property of the potential, 7^(i) can be expressed in the form 

7 f (t) = 1 [l - (-l)»W+"Ap«*(0)] , (67) 

where 7 = 7+1(0) + 7^(0), and Ap eq (0) = _p + ? i(0) — pl 9 i(0), p^{0) being the equilibrium population of the state (3 
corresponding to the rates taken at time t = 0, i.e., p e ^{0) — [J^-i 7+1(0) + fi/3,+1 7^(0)] /j. Notice that for the 
rectangular input signal in Eq. (|65|l 7 K (i) = 7 K = 7. We can also write 



A 7 K W = 7*(*) - 7 K i (*) = -(-l)" (t W 9 (0). 
As shown in Ref. [2{|, the long-time probabilities p±i(t) are given by 

1 



(68) 



P-i(t) 



!_(_!)«(*) Ap eq (0) 



,-1 \t— 



7T 

1 + e-T" 



(69) 



and p+i{t) = l-p-i(t). 

Replacing the above expressions into Eqs. (|62|) and |jf)3J), one obtains after some lengthy simplifications that 



7T7 



1 + [Ap 



eq. 



4 tanh 



IT 



(70) 



and 



= 7Tf7 out -^[A^(0)] 4 



tanh ( 

4 



^[A^(0)] 2 {l-[Ap e "(0)] 2 } 
IT 



x { 12 tanh ( ) - 7T 



1 + 2 



sech I ^-r- 



~] 2' 



(71) 



In the next subsection we will compare these analytical results for the averaged output frequency and phase diffusion 
with results obtained from a numerical solution of the stochastic differential equation Q . 



A. Comparison with numerical results 

Following the algorithm developed by Greenside and Helfand |26l Wl\ (consult also the Appendix in Ref. 28]), we 
have integrated Eq. Q for a large number of noise realizations, M, starting from one of the minima g Qo (0). From 
the initial instant of time, which we set equal to zero, we start monitoring the switches of states and recording the 
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instants of time at which those switches occur, according to Eq. {3}. We will denote by the instant of time of the 
n-th switch of state in the i-th trajectory. From the switching times i""' , the realization of N ao,a (t) corresponding to 
the i-th trajectory, AH* ' (i), can be easily calculated using Eq. (JSJ), and the corresponding realization of the discrete 
phase, (t), by Eq. 0. The noise averaged phase is then obtained by 



'(*), 



(72) 



and the phase variance by 



i M 



M 2 



A I 



$>r ,0 (*) 



(73) 



After a sufficiently long number of periods L for the system to "forget" the initial preparation, the averaged output 
frequency is calculated from the expression 



a. 



(ip a °'° [(L + 1)T]) - ((p a °<°(LT)) 



(74) 



and the averaged phase diffusion from 



D n 



v a ,0 [(l + 1)T\ 



„O!0,0 



(LT) 



T 



(75) 



Figures show the results of the numerical solution just described for a periodic rectangular input signal with 
angular frequency £1 = 2n/T — 0.01 and three different values of the driving amplitude: A — 0.14 (see Fig. 
A = 0.25 (see Fig. |3J, and A — 0.3 (see Fig. 0J. As reported previously in Ref. [17|. for high enough, but still 
sub-threshold, driving amplitudes, the present system exhibits a noise-induced frequency locking, i.e., starting from 
a non-zero value of the noise-strength D, the frequency of the output signal matches the frequency of the input 
signal, until, for strong noise, the output signal becomes desynchronized again. This effect is accompanied by a 
very pronounced suppression of the phase-diffusion of the output signal, i.e., a noise-induced phase locking. For the 
relevant values of the noise strength D, our analytical estimates Eq. ifTHjl and Eq. f7T|l agree very well with the results 
obtained from the numerical solution. Only for rather strong noise, a noticeable deviation can be observed. In this 
regime, the Kramers rates 1)561 are no longer valid. We have also plotted the results of the previous work ^tJi using 
the rates ()56t . We note that while the improvement of our analytical estimates for the frequency synchronization is 
moderate only, our prediction for the phase diffusion is strongly improved. 



VI. CONCLUSIONS 



With this work we have investigated in detail the phenomenon of frequency and phase synchronization in bistable, 
periodically driven stochastic systems. This objective is not only of foremost interest for the well known phenomena 
of Stochastic Resonance 0, and the topic of rocked Brownian motors but also carries great potential for the 
study of driven stochastic neuronal dynamics and driven excitable sytems per se |29j . Our approach takes a new look 
on this prominent problem. Starting out from a driven, Markovian continuous dynamics we derived in great detail the 
stochastic renewal dynamics of the noise-induced switching events. This contraction of the full Markovian dynamics 
in state space onto the discrete counting process of subsequent switches between the metastable states implies a non- 
Markovian dynamics for the switching times and the corresponding phase dynamics whose explicit time-evolution 
depends on initial preparation effects. The resulting non-Markovian expressions still contain the full information of 
the driven dynamics in the relevant state space and thus are not readily accessible for analytical estimates. 

In contrast, for weak noise and slow external driving, the dynamics of the underlying process simplifies consid- 
erably Consequently, in the long time limit the phase dynamics now assumes again a Markovian nature. In this 
regime we put forward new results for the phase diffusion and the frequency synchronization. In doing so, we have 
employed rectangular-shaped periodic driving signals. This choice entails two distinct advantages, namely (i) it allows 
a convenient analytic analysis of the corresponding synchronization quantities and (ii) its two-state character is also 
known to optimize the efficiency for the synchronization features. The same optimization feature holds true for related 



12 



a 



0.03 
0.02 
0.01 


-1.5 
-2 



D) 

° -2.5 



-3 











y 

X> v- 
/ 








y> x 


X 


-^^^ .I.I.I. 






yjxr' 


X 

x 


X 

X 




\ /X 



















0.02 



0.04 



0.06 



0.08 



0.1 



FIG. 2: Averaged output frequency (upper panel) and averaged phase diffusion (lower panel) as a function of the noise 
strength D for a periodic rectangular input signal [see Eq. and text below] with amplitude A = 0.14 and angular frequency 
= 2n/T — 0.01 . Solid line: Analytical results obtained from Eq. 1701 (upper panel) and Eq. 171H (lower panel), respectively. 
Dashed line: Theoretical result from Ref. [r?|. Crosses: Precise numerical results. In the upper panel, a horizontal dotted line 
indicates the frequency of the input signal. 
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FIG. 3: Like Fig.Elbut for a driving amplitude A = 0.25. 



effects such as the achievements of optimi zing the gain for Stochastic Resonance j2£| or the enhancements of energy 
transduction in driven chemical reactions |3lj| . 

Our analysis is in the spirit of prior works |l7l Il8| : our novel estimates, however, quantitatively supersede in 
accuracy those prior results, cf. the detailed comparison performed above. While the improvement for the frequency 
synchronization is moderate only, the novel estimates present a sizable improvement for the role of the phase diffusion. 
In the weak-noise regime and for slow external driving, the Markovian theory provides a very good agreement with 
numerical precise simulations. This being so, we are confident that the new insight gained into the complexity of 
non-Mar kovian, driven switching time dynamics together with its Markovian simplification obtained at weak noise 
and slow driving will prove useful for modelling and interpreting stochastic synchronization phenomena in driven 
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FIG. 4: Like Fig. but for a driving amplitude A = 0.30. 



mctastable and excitable dynamics. 
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